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We discuss the production of 771-° and 7r°7r° pairs with a large invariant mass at collider 
energies. We present a study based on a perturbative QCD calculation at full next-to-leading 
order accuracy, implemented in the computer programme DIPHOX. We give estimations for 
various observables, which concern the reducible background to the Higgs boson search in the 
channel H — > 77, in the mass range 80-140 GeV at the LHC. We critically discuss the reliability 
of these estimates due to our imperfect knowledge of fragmentation functions at high z and a 
subtle interplay between higher order corrections and realistic experimental cuts. Whereas 
the invariant mass spectrum of photon-pion pairs is theoretically better under control, in the 
dipion case large uncertainties remain. Finally we comment on the impact of our findings on 
Higgs boson searches at the LHC. We conclude that the qualitative statement that the pion 
backgrounds should not be dangerous for the H — > 77 search channel remains true at the 
next-to-leading order level. 



1 Introduction 



A major motivation for the CERN Large Hadron Collider (LHC) experiments is to shed 

light into the mechanism of electroweak symmetry breaking. As the LEP precision measurements 
favour a relatively light Higgs boson, its decay into a pair of photons is a prominent search channel 
which deserves intense consideration from the theoretical side. The inclusion of QCD radiative 
corrections is mandatory not only for the signal but also for the background if one is aiming towards 
reliable predictions. In this respect it is important to note that especially at LHC, but also at the 
Tevatron, prompt photon signals — photons which do not stem from mesonic decays — are in general 
contaminated by ir° events. The photon pairs from the mesonic decay which are very collimated, if 
the pions have a high momentum compared to the pion mass, may appear as single (fake) photons 
inside the detector. This background is especially worrisome in collider experiments, such as the 
Collider Detector at Fermilab (CDF) and DO || experiments at the Fermilab Tevatron, where the 
pions cannot all be reconstructed event by event. At the LHC jir° and 7r°7r° production represent an 
important fraction of the reducible background to the search of neutral Higgs bosons in the channel 
H — > 77 in the intermediate mass range 80-140 GeV. Since the jet-jet cross section is about eight 
orders of magnitude larger than the expected Higgs boson signal for a standard model Higgs boson, 
this reducible background is overwhelming before any selection cut is imposed. 

Regarding both issues of prompt diphoton production and Higgs boson search, various cuts, es- 
pecially isolation cuts in transverse energy have to be imposed to the candidate events, in order 
to reject most of this background. Yet, the actual fraction of this background which passes these 
cuts is not negligible, and it is not well known. Up to now its estimations rely on Monte Carlo 
simulations encoding the partonic subprocesses at the leading logarithmic (LL) accuracy only. We 
may thus worry about large higher order corrections which may be not accurately reproduced by 
parton showers. Furthermore, the hadronisation models matched with the partonic ingredients in 
these simulations are tuned to describe correctly the bulk of pion production, but they may not be 
accurate in the kinematical region relevant for the pions which pass the cut. A tentatively quanti- 
tative study of isolated jtt and 7r°7r° production, confronted to a careful analysis of data collected 
at the Tevatron could improve this situation. The motivation of this article is to provide a study of 
full next-to-leading order (NLO) accuracy, based on the computer code DIPHOX as an attempt in 
this direction. 

The computer program DIPHOX on which our study relies is a Monte Carlo code of partonic event 
generator type. This code has been designed to describe the production of pairs of particles in 
hadronic collisions, accounting for all contributing partonic processes at full NLO accuracy. These 
particles can be prompt photons or hadrons, in particular neutral pions. This code is flexible 
enough to accommodate various kinematic or calorimetric cuts. Especially, it allows to compute 
cross sections for isolated 77r° and 7r°7r° pairs, for any infrared and collinear safe isolation criterion 
which can be implemented at the partonic level. The physical content and schematical description of 
DIPHOX has beengiven in Ref . 0] , where we have concentrated on the case of diphoton production 
[|^ — n~o|l - In Ref. JllJ we recalled briefly how to use this code also in the case of 777° and 7r°7r° 
production, and we focused on the phenomenology of dipion production at fixed target energies, 
with a comparison to the data of the E706 experiment j^J at the Tevatron. A similar study was 
also presented in |l3f| . The aim of the comparison was to confront the theoretical framework on 
which DIPHOX relies to existing dipion data in order to validate this framework before transporting 
this knowledge to collider energies GJ using the factorisation property of perturbative QCD. 

The main goal of the present article is to provide a NLO prediction for the reducible background to 
Higgs boson searches in the mass range 80-140 GeV at the LHC. In section 2 we present estimates 
for distributions which are of interest regarding the Higgs boson search, such as the distributions 
of invariant mass and transverse momentum of a pair. In section 3 we will extensively discuss the 
uncertainties of our next-to-leading order approach induced by severe experimental cuts. Finally 
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section 4 gathers our conclusions and some outlook. 

2 Predictions for the 77:°, 7r°7r° background at the LHC 

In this section we discuss NLO predictions for the 771" and 7r°7r° reducible background to the neutral 
Higgs boson search channel H — > 77 in the intermediate mass range 80-140 GeV at the LHC. All 
NLO curves where computed with the DIPHOX code, version (1.2) fj]. For all figures we are using 
the MRST 2 |15|] parton distribution functions, and, if not stated differently, the KKP fragmentation 
functions [[L6| . For all plots given below the following cuts on the transverse momenta, pseudo 
rapidities of the photons/pions and their invariant mass were applied: 

Pr(7i,7Ti) > 40 GeV, 

Pt (72,^2) > 25 GeV, 

\y(l,ir)\ < 2.5, 

80 GeV < m 7r7r>77r>77 < 140 GeV. (1) 

In addition, isolation cuts are imposed to reject events which are accompanied by a certain amount 
of hadronic energy, as this is typically not the case for photonic Higgs decay. A photon is said to 
be isolated if, in a cone in rapidity and azimuthal angle about the photon direction, the amount of 
deposited hadronic transverse energy E T ad is smaller than some value Ermax fixed by the experi- 
ment: 

E* ad <E Tmax inside (y - Vl f + {<P - < R 2 (2) 

The actual experimental criteria, which are defined at the detector level, are much more complicated 
than the simple one given by Eq. (^), see p[, pfl . They cannot be implemented in a partonic 
calculation. The schematical criterion of Eq. (|2j) is a modelization of their effects at the parton 
level. The topic of isolation of photons based on the above criterion is extensively discussed in the 
literature |Q-p2|. We stress that there is a subtle interplay between higher order corrections and 
the experimental cuts. These subtleties will be discussed in the next section. 

In Fig. [l] we show the invariant mass spectra of 77, 771"° and 7r°7r° pairs for a very loose isolation 
criterion of Et max = 100 GeV in a cone of R = 0.4 at NLO. Note that the fully inclusive case would 
be plagued by uncertainties stemming from the low— z regime of the pion fragmentation functions (see 
discussion below) . To avoid this source of uncertainty we do not present figures showing the inclusive 
rate. The plot shows the importance of the reducible ir°j, ir°ir° backgrounds before isolation cuts. 
The 7r°7 component is more than an order of magnitude above the 77 background. The 7r°7r° 
background is of the order of several hundreds of pb/GeV. Tightening the isolation requirement by 
imposing Et max = 10 GeV changes the production rate of these components. This is shown in Fig. 0. 
We remind that the very narrow Higgs signal will be a sharp peak — its width is determined by the 
detector resolution of around 1 GeV — ranging from 50 to 100 fb/GeV in the given mass window pj| . 
The 7r°7r° background is reduced considerably whereas the 771° background is still of the order of the 
irreducible diphoton background. Concerning the prompt photon pairs, the used isolation criterion 
reduces its rate by about a factor of two, because the fragmentation component of the prompt photon 
pair is vetoed by isolation cuts. Note that for the given isolation criteria not only direct photon pairs 
remain. For Exmax = 10 GeV almost 20 per cent of the total rate still come from fragmentation 
in the given plot for the scale choice used. Fig. || shows the suppression of the irreducible diphoton 
background when isolation is imposed. To suppress photons from fragmentation more efficiently 

1 In this version a bug was removed compared to the first version described in In that paper the main issue was 
the prediction of photon pairs. Only fragmentation contributions in the case of severe isolation cuts were affected by 
the bug. These are numerically irrelevant for the diphoton rate predictions. 
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harder isolation cuts than the ones shown in the plot have to be imposed. A remaining uncertainty 
for the prompt photon pairs comes from higher order corrections to the box contribution. These 
corrections have been computed recently |23] but are not yet implemented in DIPHOX. To give an 
idea of the rejection factors due to isolation we have plotted the invariant mass distributions of the 
photon-pion pair for maximally allowed hadronic energies of Et max = 100, 20, and 10 GeV in a 
cone of R — 0.4 in Fig. [|. Due to the collinear fragmentation model used the hadronic energy is 
mainly located along the fragmenting particle. This leads to the fact that there is no significant R 
dependence present. Only genuine higher order 2 — > 3 matrix elements can induce such a dependence 
and we only state that they are on the per cent level. In the case of photon-pion production two 
mechanisms for the production of the photon are possible. The photon can be produced directly in a 
hard interaction or through fragmentation. In the case of isolation the latter component is reduced 
by an order of magnitude as can be seen in Fig. [5. For completeness we provide with the plot for 
the invariant mass distribution of pion pairs, Fig. 6|. In a realistic experimental situation the pion 
rates can be suppressed further. In the relevant pr range, an extra reduction factor of 2 to 3 per 
pion can be expected in fine-grained detectors with recognition of isolated tt° as two overlapping 
showers p5| , ^| . 

Another spectrum of relevance is the transverse momentum distribution of the pair produced parti- 
cles, qx = \pt i +Pt 2|- In Fig- this distribution is shown for two isolation criteria. Recall that the 
full pion photon rate is built up out of the direct component and the fragmentation component. At 
colliders which operate at lower energies the fragmentation component is typically only a fraction of 
the full rate. Through isolation the fragmentation component dies out quickly. Fig. || gives an idea 
of the suppression factors due to isolation. Note that this theoretical curve is not reliable in the first 
few bins, as multiple soft gluon emission alters the divergent behaviour of the fixed order calculation. 
Near qx = the partonic cross section diverges logarithmically. This divergence is diluted by convo- 
luting the partonic cross section with the pion fragmentation function. Also at qx = Et max there is 
a critical point inside the physical spectrum (see for a discussion). The large bin width prevents 
from observing this singular behaviour. To deal with these issues in a rigorous way is beyond the 
scope of this study0. 

Hard isolation cuts test fragmentation functions at high z values, where current parametrisations 
are not at all constrained by experimental data. This problem is visualised in Fig. ^| where we 
show the photon-pion invariant mass distribution calculated with both, the KKP and the older 
BKK |27]] fragmentation functions. The KKP parametrisation leads to a higher prediction than 
the BKK one. The ratio of the two curves (KKP/BKK) is about 3/2. For dipion observables the 
situation is much worse. In Fig. |l0| the discrepancy between the BKK and KKP parametrisations in 
the critical region is shown. One observes that at high z values the fragmentation gluon-to-pion is 
strongly suppressed. Focusing on the quark-to-pion fragmentation, one sees that for high z values 
the parametrisations differ considerably. This will be discussed in more detail below. 

We have to comment now on the remaining theoretical uncertainty due to higher order corrections. 
As will be pointed out below the large z region is dangerous for the predictive power of NLO 
calculations because terms are produced which are divergent as z — ► 1. In Fig. |ll| , we plot LO and 
NLO predictions^] for the invariant mass distribution of the photon-pion pair for several isolation 
criteria. Note that the ratio of both curves is pretty stable with decreasing Et max ■ The ratio of 
both curves lies around 1.7, a typical value for a NLO correction. We just mention that the situation 
worsens as Et max 

is reduced further which indicates the presence of large NLO contributions induced 
by pressing z to one. As will be pointed out below there are two reasons for that. Phase space regions 
which are excluded at LO are populated at NLO and become dominant, and large logarithms of 

2 The corresponding resummations to be done are technically very involved and no satisfying procedure is defined 
yet. 

3 LO refers here to the Born matrix elements whereas NLO parton densities and fragmentation functions are used. 



3 



(1 — z) are introduced. To have a proper control on the latter fixed order effect, a resummation of 
the relevant contributions is needed. For the isolation criteria used in the given plots we conclude 
from the NLO to LO ratios that these effects are not yet dominant. Still, for harder isolation criteria 
the predictivity of our fixed order approach has to be critically investigated. 

As the isolation cuts spoil compensations between scale dependent LO and NLO contributions 
the scale dependence remains important. To give an idea of the uncertainties we have plotted 
the variation of the invariant mass spectrum of 7r°7 pairs under the change of renormalization 
(/i), factorisation (M) and fragmentation scale (Mf) for a strict isolation criterion in Fig. [l2]. By 
multiplying our standard scale choice by 1/2 and 2 one finds variation of about +40 per cent to -30 
per cent. This variation is reduced if one uses less restrictive experimental cuts. 

The importance of NLO calculations is best visible in distributions which are restricted by LO 
kinematics. No LO parton shower can provide for a correct description of events far away from the 
LO kinematics. A relevant observable to exemplify this point is the qr distribution of a photon-pion 
pair. Increasing the isolation criterion, suppresses the LO phase space more and more, such that 
finally the whole tail of the distribution is a pure NLO effect, see Fig. |l3]. For completeness we show 
the effect of isolation on this observable in Fig. For a comparison of the shape of our DIPHOX 
prediction to PYTHIA see §§. 



3 Experimental cuts and higher order corrections 

To understand the reliability of the presented predictions it is necessary to discuss the interplay 
between experimental cuts and theoretical NLO results. QCD cross sections are stabilised by the 
inclusion of higher order corrections, as scale dependencies are reduced. This works generally the 
better the more inclusive the observables are. Yet, the collider experiments at the Tevatron and 
the forthcoming LHC do not perform inclusive photon measurements. The experimental selection 
of prompt photons requires isolation cuts. We always use the cone isolation criterion introduced in 
Eq. (||) . For the hadronic energy propagating into the isolation cone two sources can be distinguished 
in an experimental setup. First there is a component which stems from the actual partonic, hard 
reaction which may be mapped to certain Feynman diagrams - this component is calculable inside 
a QCD calculation to a given order. And second there is a component coming from multiple- 
interactions and other events in the same particle bunch (pile-up). One has to be aware of the fact 
that an experimentally used Exmax is thus always bigger than a partonic value used in our study. 
It is even so that the experimental values Et max are so small that a large fraction of the allowed 
hadronic transverse energy accompanying the candidate photon may come from such events. This 
leaves not much phase space for partonic fragmentation reactions. In the following we mean with 
Ermax only the hadronic energy stemming from the hard interaction. The sensitivity on Ermax 
defined in Eq. (|2|) is best seen by looking at the fragmentation variable z, which is the momentum 
fraction of a hadron, ptk (h = tt,7), and its parent parton pt- At LO the kinematical constraints 
on this variable are induced by the experimental cuts through the collinearity assumption in the 
fragmentation model, pr — zpr + (1 — z)pr = prh + E^ ad . One finds 

_ PTmin 

" ^ ^min — i TP V / 

PT min + HjT max 

where prmin is the experimental cut on pru- Evidently, the variable z is pressed towards one, as 
the ratio Et max / 'pt min decreases, means for hard isolation. At NLO, where extra radiated partons 
are present, this bound is even more severe as the veto is more easily fulfilled. The bulk of the pions 
which are produced in collider events correspond to values of z substantially smaller than 1. Hard 
isolation tests only the very end of the z spectrum. The question is, wether we can reliably estimate 
the production rate of these high-z pions, and especially of hard jir° and ir°ir° pairs. We find three 
aspects important in this respect, which are discussed in more detail in the following. 
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3.1 Reliability of fragmentation functions 

The usual fixed order perturbative QCD calculations of hadro production - as e.g. the one coded in 
DIPHOX - rely on a model of independent collinear fragmentation through fragmentation functions 
of a parton of species j (quark of any flavour, or gluon) into a pion, Dj^^o (z, Mj), z is defined 
above. Let us first comment on the uncertainties at small z, where this treatment is questionable. In 
particular the fixed order calculations do not resum the In (1/z) arising order by order in perturbation 
theory, which become large in the small z region. The latter are associated with coherence effects in 
the multiple emission of soft gluons, which control the inclusive production of small z pions ]29|,|3Cj|]. 
Therefore this usual framework is notoriously inadapted to account for the production of these 
small z pions. Moreover the dominant z values at colliders are smaller than for fixed targets, as a 
consequence of scaling violations. Therefore the erroneous treatment of the fragmentation process 
at small z will plague the estimates of inclusive dipion observables at colliders. This is the reason 
why the fits of independent fragmentation functions D^^o (z, M"j) to experimental data presented 
in the literature |jqp7|, and which we use in this work, start above some minimum value for z, 
e.g. z m i n ~ 0.1. The independent fragmentation framework is expected to account correctly for the 
production of pions with larger z values. The requirement of a loose isolation criterion such as, for 
example Exmax < 100 GeV for pions with pt-k > 25 GeV imposes that z > 0.2 over the whole 
Pttt range of interest. To avoid uncertainties due to small z values, as was pointed out already, we 
do not provide estimates for inclusive observables, but only for observables corresponding to pions 
submitted to isolation cuts. 

The behaviour near the end point, z — ► 1, is responsible for most of the uncertainties that plague 
the present study. Indeed, the use of very tight isolation cuts such as those used by the collider 
experiments select the vicinity of the elastic boundary z ~ 1. For example, in the Run lb of the 
Tevatron, for prompt photon measurements, CDF required Exmax — 2 GeV in a cone R — 0.7. 
For pr > 15 GeV, this means z m i n (pT) > 0.9. The value Ermax = 2 GeV required by the CDF 
experiment includes also the deposition of hadronic transverse energy coming from underlying events, 
while the latter are ignored in our partonic calculation as noted above. Consequently, the actual 
effective value of E^max which should be used in a partonic calculation should be even smaller, the 
effective z m i n would be even larger. However the fragmentation functions Dj_> n o (z, M?) are fitted to 
data corresponding essentially to z < 0.7, i.e. much below the relevant region for the present study. 
Moreover, the data used in jl6|,[27|] concern hadronic final states in e + e~ annihilation. The gluon-to- 
pion fragmentation function is poorly constrained at moderate and large z by this process because 
the subprocesses involving gluons in e + e~ annihilation appear at higher order in the perturbative 
expansion in powers of a s , while subprocesses involving outgoing gluons contribute at lowest order 
in hadronic collisions. Although the contributions from gluon fragmentation are not expected to 
dominate, large uncertainties on their magnitude may still affect quantitative predictions. In order 
to illustrate how poor the present-day control over the large z region is, it is instructive to compare 
the extrapolations of the KKP and BKK sets of parton-to-7r fragmentation functions. The KKP set 
is the most recent one, and it is constrained by more data, especially the LEP1 data at the Z peak, 
in the region 0.1 < z < 0.7. These two sets use different analytic ansatze. These different ansatze do 
not differ very much in the z range in which they are most constrained by the e + e~ data of course, 
but their respective extrapolations turn out to differ substantially from each other above z > 0.7, 



especially between z = 0.9 and 1, as shown in Fig. 10. This is true for the gluon, but also, though to 



a lesser extend, for the u and d fragmentation functions. This translates into a large uncertainty for 
the predictions of the reducible Higgs boson background in the channel H — > 77 in the mass range 
80-140 GeV. Clearly a determination of fragmentation above z > 0.7 from data sensitive to this 
range is required for any further accurate prediction. Such an experimental input could in principle 
be provided by data on inclusive ir° production at fixed target energy |Il"|fi as well as measurements 

4 In Rcf. |3l| , the comparison between E7 ri data and BKK fragmentation functions reveals an overall normalisation 
discrepancy. However, as noticed in Ref. [JLl|, this problem seems to be largely reduced by using the KKP set. 
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of isolated pions by the collider experiments CDF and DO, in association with their prompt photon 
measurements. 



It is worth mentioning that the present inaccuracy of the description of fragmentation at large z is 
not only a problem confined to the type of perturbative QCD calculations which are discussed in this 
article. They also potentially affect the fragmentation stage of the "full hadronic event generators" 
such as PYTHIA 3^| and HERWIG J|3| as well. Indeed, the respective fragmentation ingredients 
of these event generators have been tuned to describe accurately the bulk of the data, especially at 
LEP, whereas the tails of the distributions near the end point have, to our knowledge, never been 
analysed in detail. In addition, the parton showers implemented in these event generators provide 
an effective summation of the soft gluon effects near z ~ 1 at a leading logarithmic accuracy; on the 
other hand they miss a part of the higher order corrections to the partonic hard subprocesses which 
may be large. The respective approximations in the "partonic calculations" on one hand, the full 
event generators on the other hand, being different, it would be instructive to compare the results 
given by both approaches pgfl. 



3.2 Large z values and the need for end-point summation 

As the production of isolated pions at colliders and the one of inclusive pions at fixed target ex- 
periment, deal with z values near the end point z = 1, the theoretical calculation of the partonic 
subprocesses is confronted with the appearance of powers of large logarithmic terms In (1 — z) to 
any fixed order in perturbation theory. These large logarithms have in principle to be summed to 
all orders in order to improve - if not restore - the predictive power of the partonic calculation. The 
general framework to resum infrared logarithms in hadronic collisions has been settled in |54],|35|], 
however further work is needed to treat the endpoint summation in the case of one particle pro- 
duction, even more for pair production, where the kinematics is more constrained. In the present 
study, we investigate what can be done without such summations, relying on the extrapolations of 
the fragmentation functions available in the literature, and examining critically the limitations of 
our approach. 

To be explicit let us consider the process PP — ► jir a + X in detail. The photon can be either direct 
or from fragmentation. The second case is considerably suppressed if one imposes severe isolation; 
thus we deal with the direct process only. At leading order there arc two types of partonic processes: 
qg — ► qq — + 717. At small x values typical for the LHC the first is dominant. Radiative 

corrections give rise to initial state and final state singularities. Only the latter are of relevance for 
the point to be made here. They arise from the splitting of the outgoing parent quark into a quark 
and a gluon. Both may fragment into a pion but, again, for severe isolation one may focus on the 
dominant q — > 7r° fragmentation only. Thus we are looking at the case where a collinear gluon is 
radiated of a hard quark, indicated by the subscript g / / q, which subsequently fragments into a pion. 

5 Here q stands for quark or antiquark. 
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Following Ref. H the finite remnant of the collinear approximation of this correction is of the form 

cr(PP-> qgi)\ g //q = 



dyndy-fdpTiT / — 



1 1 

dx f dz a s (n) 



z 2?r 



xD^o /q {x,M 2 f 
'ln(l - z) 



F q/P { Xl ,M 2 ) F g/P (x 2 ,M 2 ) 




(4) 



4%) 



4? 0) 



(1 



+ -C F S(1 - z) 



C F {l + z 2 ) 



Here y*,pr n are the rapidity and transverse momentum of the pion, y y the rapidity of the photon, 
xz the momentum fraction of the pion with respect to the parent quark, pt-k — xzpxq, and z 
the momentum fraction of the gluon with respect to the parent quark, Pt g — zpx q - In the given 
collinear approximation the momentum fractions of the initial partons are 



■C2 



Pt , 



Pt 



Vs 

PTq+PTg 



~[exp(y n ) +exp(y 7 )] 
[exp(-y^) +exp(-y 7 )] 



The value of x m i n will be discussed below, as it is the sensitive variable to isolation cuts. The 
dots stand for terms which are less dominant as z — > 1. They are not necessarily small. By writing 
aqq\z)/z = Cf[(1 — z) 2 /z+2], one sees that the first term in the bracket only leads to non logarithmic 
corrections. Thus, approximating the bracket by 2 and using that for z — > 1, ln(z) = 0(1 — z), the 
z integration can be done explicitly and one finds 



a{PP -> qgi)\ g // q 



dy^dy^dpr 7T 



dx a s (p)Cf 



x£) 7r o / ,(x,M|) 



F q/P ( Xl ,M 2 )F g/P (x 2 ,M 2 ) 



ln(l 



Xi 

2 ' In 



X2 



\M 



LO 

qg^q-y\ 



x 2 M 2 



ln(l 



(5) 



The important point is that the behaviour of these terms changes drastically as one goes from 
the inclusive case, where x > x m i n — 2pt-k/VS is much smaller than one, to the case of severe 
isolation which press x to values near one. Using the fact that the fragmentation function behaves 
as D^o / q (x) ~ N(l — x) 13 with (3 ~ 1, one can write 



in i a ) 



l - 



+ 0(1 - x min , 1 - 0) 



(6) 
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and perform the integral over x. Focusing on the double logarithmic terms, one finds 



a{PP^ qgl )\ 



g/h 



dy^dy-dpr „ 



x(l 
1 



D^ /q {x min ,Mj) F q/P ( Xl ,M 2 ) F g/P (x 2 ,M 2 ) 



■i'2 



ln(l 



2 + ln 



ln(l 



(7) 



the dots stand for single logarithmic terms and terms which are suppressed by an extra factor of 
(1 — x min ) and (1 — f3). Note that outside the bracket, up to the factor a s (fi)Cp /(2tt), one has just 
the approximation for the leading order contribution. This means that the bracket itself times this 
factor could be taken as a good approximation for the higher order part of a if factor^. 



K = 1 



ln(l 



2 + ln 



Ptjl 



ln(l 



+ 



(8) 



In here a factor of two is included from the second process with initial partons qg. To get a feeling 
for the typical scales in the case of isolation one may try to choose Mf such that the large logarithms 
compensate for a pm ~ PTmin- This can be achieved by a (optimal) scale choice 



Mj Pt — PTn rain \ 



Et, 



PTi 



(9) 



Plugging in p Tnu 
to choose for Mf 



n =25 GeV[]andi; Tm( 

in the figures presented not the canonical choice M^ 7 /2 but rather M w 



10 GeV one finds M° f pt 



10 GeV. This motivated us 
Note 



that for a much higher scale choice the logarithms become constructive. For the upper choice of 
parameters and using Mf = 2pTmin arbitrarily large K factors are induced in the limit x m i n —> 1- 
The fact that the optimal scale in Eq. (^|) goes to zero, as Exmax goes to zero indicates that the 
cone isolation criterion is not infrared safe in this limit. In other words, if the experimental value 
of Et max is already saturated by events from the hadronic environment such that there is no phase 
space left for the hard parton, fixed order calculations are not applicable. 

It should be stressed that the whole discussion was neglecting constant terms and sub-leading loga- 
rithms. Thus it should only give a rough idea of the much more involved, complete situation. Using 
the same approximations the exercise can be worked through keeping also non-leading terms. The 
result is that one encounters in realistic situations additive contributions which are of the same size 
as the considered ones. All this deserves further investigations especially in the light of a possible 
resummation of the large logarithms and the role of these terms. 



3.3 Hard isolation and asymmetric cuts 

In LHC studies concerning Higgs boson search the used px cuts are commonly asymmetric. We want 
to make the point here that this jeopardizes the predictive power of NLO calculations. The interplay 
between such asymmetric cuts and hard isolation criteria leads at LO to vetoed phase space regions 
which are filled by NLO processes only. The effect is visualised in Fig. |lj, where the LO and NLO 

6 Note the way this K factor depends on the experimental cuts applied and on the scale choice. This is a nice 
example of the fact that for non inclusive quantities the naive use of K factors is not at all adequate to take into 
account higher order corrections. 

7 Due to the NLO kinematics Pt~/ > Pt-k- F° r asymmetric px cuts the photon will always be the particle restricted 
by the larger px cut. 
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DIPHOX prediction for the dipion pt spectrum is plotted for symmetric, both pr > 25 GeV, and 
asymmetric cuts, pxi > 40 GeV, px2 > 25 GeV. In the case of asymmetric cuts one observes that 
the LO component is suppressed for low px values. This can be understood as follows: 
Due to transverse momentum conservation, one has 

Pt i + E T \ d =p T2 + E%tf (10) 

where the pion 1 is assumed to have higher px than pion 2. This means for the transverse momentum 
of the pair, qr that 

< <7t = Pti ~ PT2 — E T2 — E T1 < E T (11) 

and because pti > Pti min (= 40 GeV in our case) 

p T2 >p T1 -EJp a *>p Tlrnin -EJp a * (12) 

This means that in the case of px i min — E™ ax > PT2 min one has a phase space restriction coming 
from the isolation criterion. Note that a saturation of the bound, PT2 — Pti min — E T lax is only 
possible if E T °f — + 0, i.e. z\ — ► 1. But in this limit the fragmentation function goes to zero and 
suppresses the cross section additionally. This will not happen in the case of photon-pion production 
where the same formulas apply with E^°f — 0. Eq. (|lj) is also true and leads to a restriction of the 
phase space but there is no additional suppression. The problem aggravates for lower ETmax values 
but becomes less relevant when ETmax > PTmini — PTmin2 and vanishes in the case of symmetric 
cuts PTmini = PTmin2- The hole in the pt spectrum will be filled with events from higher order 
corrections. As the low pr region is energetically preferred, one gets numerically important NLO 
contribution to the cross section in the dipion case. The ratio between LO and NLO can become 
very big due to this effect. As pointed out, in photon-pion production the LO order contribution is 
less suppressed and the effect is milder. In the qr distribution shown in Fig. ||, one observes that 
only the first bin contains the LO contribution due to the same kinematical restriction mentioned 
above. On the other hand the NLO contribution is basically zero there (even below the plotted 
range). This is because pions with pt's as large as the photon pt are disfavoured. The steep fall of 
the quark-to-pion fragmentation functions makes it very unlikely to find a photon-pion pair with 
PT-y ~ PTix > 40 GeV which is needed to have ~ 0. 

From the discussion one sees that symmetric cuts are preferable from theoretical side for processes 
which have to be suppressed by isolation. Whereas this effect is for sure important for dipion 
production the effect is numerically sizable in the photon-pion case only for isolation criteria which 
are harder than the ones given in the presented figures. As can be inferred from Fig. |ll|, the 
ratio of the NLO to LO curves does not explode because the discussed phase space effect is not 
yet numerically sizable. Still we want to make the point that our NLO analysis suggest to choose 
symmetric transverse momentum cuts as theoretical higher order predictions tend to be more reliable 
in that case. 

4 Conclusions and outlook 

We have presented a full next to leading order study of the production of 7r°7 as well as tt° pairs 
with a large invariant mass relying on the computer code DIPHOX. We have presented estimates 
for invariant mass and transverse momentum of the pair {qr) distributions taking into account 
realistic experimental cuts. The NLO corrections are large. Theoretical calculation improves the 
scale dependences in the case of inclusive observables. On the other hand, the results remain plagued 
by large uncertainties even in the NLO approximation, if realistic isolation criteria are applied. One 
source of large uncertainties is our poor present knowledge of the parton-to-pion fragmentation 
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functions at large z. The estimates which we have presented here rely on extrapolations of analytical 
ansatze for these fragmentation functions for values of z above 0.8, i.e. outside the region where they 
were fitted to the data. The spreading of the extrapolations of different sets of parametrisations 
available in the literature for z above 0.8 illustrates this fact. Experimental information may be 
provided by inclusive measurements of single pions and pions pairs at fixed targets, e.g. as performed 
by the E706 experiment [12], as well as by measurements of isolated pions at colliders, as could be 
performed by the CDF ]^| and DO || experiments. 

Apart from this more experimental issue we pointed out that our next-to-leading order approach is 
still plagued by relatively large uncertainties due to phase space effects related to the different LO 
and NLO kinematics and also the presence of potentially large logarithms ln(l — x m i n ), if severe 
isolation cuts are imposed. Concerning the phase space effect, we discussed that symmetric pr cuts 
for the pions/photons are preferable from theoretical side. Otherwise NLO corrections are magnified 
in phase space regions vetoed at lowest order. Our predictions for 777° production seems still reliable, 
up to scale variations of around 40 per cent, for the used isolation cuts. The same cannot be said for 
7r°7r° production. The region z —* 1 is probed in a more critical way in this case. Further, predictions 
are very sensitive to the unconstrained part of fragmentation functions. We have to conclude that 
the uncertainties for this component of the reducible background for Higgs searches are not under 
control from theoretical side even at the next-to-leading order level. We only quote that we find 
large K factors and large scale variations, factors of 3 to 4 respectively, if we use the same isolation 
cuts as was done for the isolated 771- production case. Fortunately this component to the diphoton 
background can be suppressed experimentally by separating pions from photons by distinguishing 
the shapes of the corresponding electromagnetic showers. Reduction factors of 2 to 3 per pion are 
possible and compensate somewhat the poor knowledge of the pion rates. Still, a more satisfactory 
theoretical knowledge of the dipion and pion/photon rates at LHC requires a better understanding 
of fragmentation functions near z = 1, the resummation of large endpoint logarithms of (1 — z), 
together with the use of experimental cuts which do not destabilise compensations between LO and 
higher order contributions. 

From the discussion above it should be clear that it is a delicate issue to make definite predictions 
for the 77T° and 7r°7r° component of the H — ► 77 background. All statements depend strongly on 
the experimental cuts. Especially severe isolation cuts restrict the predictivity of our result. To 
make some explicit statements we just want to remind how the 77r , 7r°7r° rates compare to the 
irreducible 77 background for a definite isolation criterion, see Fig. |[ As the Higgs boson of the 
Standard Model is already restricted by the LEP experiments to have a mass of Mh > 115 GeV, 
the 77T° rate is about a factor of two below the irreducible 77 background. Uncertainties due to 
fragmentation functions and scale variations may be factors of around two. On the other hand 
additional experimental capabilities based e.g. on calorimetric studies indicate that the rejection 
factor of pions w.r.t isolated photons can be improved by at least a factor two [[j],[^]. These numbers 
suggest a conservative upper bound on the ratio of the 771-°, 77 rates of 0.5. If one tightens the 
isolation criterion this will be further reduced, though one has to keep in mind that the theoretical 
uncertainty is getting more important. In this respect we note that a comparison with PYTHIA 
distributions shows that our rates are lower if compared to default PYTHIA settings |28| . Though 
we want to stress the relevance of a next-to-leading order analysis, from our discussion one may 
not expect that the 771"° process becomes a dangerous background. As already noted above the 
uncertainties from fragmentation functions and scale variation for the 7r°7r° rate are substantial. To 
be conservative factors as big as ten should be considered. Assuming again a reduction factor of 
two per pion from pion-photon separation would lead in Fig. ^| to a 7r°7r° to 77 ratio of around 
1/3 for M-xti > 115 GeV. More severe isolation will suppress this ratio considerably again with the 
price of even increased uncertainties. Note that the given numbers should only give an idea of the 
situation, as they depend strongly on the experimental settings. We can tentatively conclude that 
next-to-leading order corrections should not change the qualitative statement that the 7r°7r°, 7r°7 
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backgrounds are sufficiently under control at the LHC, although a couple of interesting QCD issues 
plaguing them still deserve to be addressed. 
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Figure 1: Comparison of the 77, 771- and ir°ir° invariant mass distribution for the isolation cuts 
E T ma X = 100 GeV, R=0.l 
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Figure 2: Comparison of the 77, 771- and 7r°7r° invariant mass distribution for the isolation cuts 
E Tmax = 10 GeV, R=0.l 
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Figure 4: Invariant mass distribution of photon-pion pairs for various values of Ermax- 
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Figure 5: The suppression of the fragmentation component of tt°j production through isolation. 
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Figure 7: The qr dependence of photon-pion pairs split into direct and fragmentation part with (top) 
loose isolation and (bottom) severe isolation. 
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Figure 8: The qx dependence of photon-pion pairs for various isolation criteria. 
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Figure 9: Comparison of the KKP and BKK fragmentation functions for ir° 7 production. 
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Figure 11: Comparison of LO and NLO predictions of the invariant mass distribution of photon-pion 
pairs for various isolation criteria. 
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Figure 12: Scale dependence of the incvariant mass distribution of photon-pion pairs for the isolation 
criterion ETmax = 10 GeV in the cone R = 0.4. The scales are varied around our standard choice 
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Figure 13: LO and NLO contributions for the qr distribution in 7r°7 production for different isolation 
criteria. The LO and NLO contributions populate different phase space regions when severe isolation 
is applied. 
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Figure 14: The transverse momentum spectrum of the pions in dipion production for symmetric 
(bottom) and asymmetric (top) pt cuts. Scales: n = M = M 7r7 /2, Mf — M nl /8. 
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